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Assistant Commissioner for Patents October 24, 2C 

Washington, DC 20231 

Sir: 

The following preliminary amendments and remarks are respectfully submitted in 
connection with the above-identified application. 
IN THE ABSTRACT OF THE DISCLOSURE 

Please replace the original Abstract with the attached revised Abstract. 
IN THE SPECIFICATION 

Please replace the Specification with the attached Substitute Specification attached 

hereto. 

IN THE CLAIMS 

Please amend the claims as follows: 

1 . (Amended) A method for designing a technical system, comprising: 



providing a target function of the technical system, influenced by a parameter 



vector including a plurality of parameters; 
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providing secondary conditions of the technical system to determine a domain 
for valid working points of the technical system; 

determining a Stratonovich equation using a projection operator; 

solving the Stratonovich equation numerically, and thereby determining an 
efficient working point; and 

using the determined efficient working point to design the technical system. 

2. (Amended) The method as claimed in claim 1, wherein the Stratonovich 
equation is determined such that its solution varies in a domain of valid working points. 

3. (Amended) The method as claimed in claim 1 ? wherein the Stratonovich 
equation is solved numerically by a predictor-corrector method. 

4. (Amended) The method as claimed in claim 3 , wherein the predictor is 
determined with the aid of a trapezium rule. 

5. (Amended) The method as claimed in claim 3, wherein the corrector is 
determined by a Newton method. 

6. (Amended) The method as claimed in claim 1, further comprising: 
additionally performing step size control. 

7. (Amended) The method as claimed in claim 6, wherein the performance of the 
step size control includes at least one of, 

carrying out a step size reduction if a sufficiently large approximation error is 
determined, and 

carrying out a step size enlargement if a sufficiently small approximation error 

is determined. 
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8. (Amended) The method as claimed in claim 1, wherein the Stratonovich 

equation has the following form: 

t 

X t = X a - Jp(X s ) • F(X S ) • ds + Z , wherein t > a and X a = x 0j 

a 

and wherein, 

a denotes a starting instant (a > 0), 

xo denotes a permissible starting parameter vector, 

X t denotes a solution parameter vector at the instant t, 

P(X S ) denotes a projection matrix in the parameter vector X s , 

F(X S ) denotes an ascending direction, and 

Z denotes a random variable. 

9. (Amended) The method as claimed in claim 8, wherein the ascending 
direction is determined by: 

F= Vf(X s ), 

wherein, f denotes the target function. 

10. (Amended) The method as claimed in claim 8, wherein the random variable is 
determined by 

t 

Z=s ■ Jp(X s ) o dB s> 

a 

wherein, 

s denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian movement. 
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11. (Amended) The method as claimed in claim 1, further comprising: 

using the determined working point as a starting point for a local optimization 

method. 

12. (Amended) The method as claimed in claim 1, further comprising: 
transforming a secondary condition which is present in the form of an 

inequality, into an equation using a slack variable s. 

13. (Amended) The method as claimed in claim 12, wherein the secondary 
condition hj(x), which is present in the form of the inequality hi(x)<0, is formulated as 

hi(x, s) = hi(x) + g(s), 

wherein, 

g(s) denotes a function including a value range of all real values greater 
than or equal to zero, 

s denotes the slack variable, and 

hj(x) denotes a secondary condition hi(x) < 0, 

14. (Amended) The method as claimed in claim 13, wherein the function g(s) is 
determined as 

g(s) = s 2 . 

15. (Amended) The method as claimed in claim 1, wherein the design is at least 
one of a new design, an adaptation and a control of the technical system. 

16. (Amended) An arrangement for designing a technical system, comprising: 

a processor unit set up to perform the various steps wherein, the technical 
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system includes a target function which can be influenced by a parameter vector including a 
plurality of parameters and secondary conditions of the technical system are used to 
determine a domain for valid working points of the technical system, the processor unit being 
set up to, 

determine a Stratonovich equation using a projection operator; 
solve the Stratonovich equation numerically and thereby determine an 
efficient working point; and 

use the determined efficient working point to design the technical 

system. 

Please add the following new claims: 

17. The method of claim 1, wherein the secondary conditions are associated with 
the target function. 

18. The arrangement of claim 16, wherein the secondary conditions are associated 
with the target function. 

1 9. The method as claimed in claim 2, wherein the Stratonovich equation is solved 
numerically by a predictor-corrector method. 

20. The method as claimed in claim 4, wherein the corrector is determined by a 
Newton method. 

21. The method as claimed in claim 2, further comprising: 
additionally performing step size control. 

22. The method as claimed in claim 21, wherein the performance of the step size 
control includes at least one of, 
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carrying out a step size reduction if a sufficiently large approximation error is 
determined, and 

carrying out a step size enlargement if a sufficiently small approximation error 

is determined. 

23. The method as claimed in claim 2, wherein the Stratonovich equation has the 
following form: 

t 

X t = X a - Jp(X s ) • F(X S ) * ds + Z , wherein t > a and X a = x 0j 

a 

and wherein, 

a denotes a starting instant (a > 0), 

x 0 denotes a permissible starting parameter vector, 

X t denotes a solution parameter vector at the instant t, 

P(X S ) denotes a projection matrix in the parameter vector X s , 

F(X S ) denotes an ascending direction, and 

Z denotes a random variable. 

24. The method as claimed in claim 23, wherein the ascending direction is 
determined by; 

F= Vf(X s ), 

wherein, f denotes the target function. 

25. The method as claimed in claim 9, wherein the random variable is determined 

by 

t 

Z = s • Jp(X s ) o dB s , 
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wherein, 

s denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian movement. 

26. The method as claimed in claim 23, wherein the random variable is 
determined by 

t 

Z = 8 • Jp(X 8 ) o dB s , 

a 

wherein, 

c denotes a scaling factor, and 

{B $ : s > a} denotes an n-dimensional Brownian movement. 

27. The method as claimed in claim 24, wherein the random variable is 
determined by 

t 

Z = s • Jp(X s ) o dB s , 

a 

wherein, 

s denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian movement. 

28. The method as claimed in claim 2, further comprising: 

using the determined working point as a starting point for a local optimization 

method. 

29. The method as claimed in claim 2, further comprising: 

transforming a secondary condition which is present in the form of an 
inequality, into an equation using a slack variable s. 
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30. The method as claimed in claim 29, wherein the secondary condition hi(x), 
which is present in the form of the inequality hi(x)<0, is formulated as 

hi(x, s) = hi(x) + g(s), 

wherein, 

g(s) denotes a function including a value range of all real values greater 
than or equal to zero, 

s denotes the slack variable, and 

hi(x) denotes a secondary condition hj(x) < 0. 

3 1 . The method as claimed in claim 30, wherein the function g(s) is determined as 
g(s) = s 2 . 

32. The method as claimed in claim 2, wherein the design is at least one of a new 
design, an adaptation and a control of the technical system. 

33. The arrangement of claim 16, wherein the processor unit includes, 
a central processor; 

a memory; and 

an input and output unit. 

34. The arrangement as claimed in claim 16, wherein the Stratonovich equation is 
determined such that its solution varies in a domain of valid working points. 

35. The arrangement as claimed in claim 16, wherein the Stratonovich equation is 
solved numerically by a predictor-corrector method. 

36. The arrangement as claimed in claim 35, wherein the predictor is determined 
with the aid of a trapezium rule. 
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37. The arrangement as claimed in claim 16, wherein the Stratonovich equation 
has the following form: 

t 

x t = x a - J p ( x s> • p ( x s) • ds + z , wherein t > a and X a = x 0 , 

a 

and wherein, 

a denotes a starting instant (a > 0), 

x o denotes a permissible starting parameter vector, 

denotes a solution parameter vector at the instant t, 
P(X S ) denotes a projection matrix in the parameter vector X s , 
F(X S ) denotes an ascending direction, and 
Z denotes a random variable. 

38. The arrangement as claimed in claim 37, wherein the ascending direction is 
determined by: 

F= Vf(X s ), 

wherein, f denotes the target function. 

39. The arrangement as claimed in claim 37, wherein the random variable is 
determined by 

t 

Z = s • j*P(X B ) o dB s , 

a 

wherein, 

8 denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian movement. 
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40. The arrangement as claimed in claim 1 6, wherein the design is at least one of a 
new design, an adaptation and a control of the technical system. 

41. The method of claim 1, wherein the projection operator takes account of the 
secondary conditions. 

42. The method of claim 1 7, wherein the projection operator takes account of the 
secondary conditions. 

43. The arrangement of claim 16, wherein the projection operator takes account of 
the secondary conditions. 

44. The arrangement of claim 1 8, wherein the projection operator takes account of 
the secondary conditions. 

45. The method of claim 1, wherein the efficient working point is an optimum 
working point. 

46. The arrangement of claim 1, wherein the efficient working point is an 
optimum working point. — 

REMARKS 

Claims 1-46 are now present in this application, with new claims 17-46 being added 
by the present Preliminary Amendment. It should be noted that the amendments to original 
claims 1-16 of the present application are non-narrowing amendments, made solely to place 
the claims in proper form for U.S. practice and not to overcome any prior art or for any other 
statutory considerations. For example, amendments have been made to broaden the claims; 
remove reference numerals in the claims; remove the European phrase "characterized in that"; 
remove multiple dependencies in the claims; and to place claims in a more recognizable U.S. 
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form, including the use of the transitional phrase "comprising" as well as the phrase 
"wherein". Other such non-narrowing amendments include changing the phrase "or" to --at 
least one of--, and reorganizing method (using separate clauses beginning with "-ing" verbs) 
and apparatus-type claims (setting forth elements in separate paragraphs and re-ordering 
elements) in a more recognizable U.S. form. Again, all amendments are non-narrowing and 
have been made solely to place the claims in proper form for U.S. practice and not to 
overcome any prior art or for any other statutory considerations. 

SUBSTITUTE SPECIFICATION 

In accordance with 37 C.F.R. §1.125, a substitute specification has been included in 
lieu of substitute paragraphs in connection with the present Preliminary Amendment. The 
substitute specification is submitted in clean form, attached hereto, and is accompanied by a 
marked-up version showing the changes made to the original specification. The changes have 
been made in an effort to place the specification in better form for U.S. practice. No new 
matter has been added by these changes to the specification. Further, the substitute 
specification includes paragraph numbers to facilitate amendment practice as requested by the 
U.S. Patent and Trademark Office. 

CONCLUSION 

Accordingly, in view of the above amendments and remarks, an early indication of the 
allowability of each of claims 1-46 in connection with the present application is earnestly 
solicited. 
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Should there be any outstanding matters that need to be resolved in the present 
application, the Examiner is respectfully requested to contact Donald J. Daley at the 
telephone number of the undersigned below. 

If necessary, the Commissioner is hereby authorized in this, concurrent, and future 
replies, to charge payment or credit any overpayment to Deposit Account No. 08-0750 for any 
additional fees required under 37 C.F.R. § 1.16 or under 37 C.F.R. § 1.17; particularly, 
extension of time fees. 



Respectfully submitted, 



HARNESS, DICKEY & PIERCE, P.L.C 



By: 



Donald J. Daley, Kdg. No. 34,313 




DJD:kna 



P.O. Box 8910 
Reston, Virginia 20195 
(703) 390-3030 
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ABSTRACT OF THE DISCLOSURE 

A method for designing a technical system is specified in the case of which the 
technical system includes a target function which is influenced by a parameter vector having n 
parameters. Furthermore, secondary conditions are prescribed with the aid of which a domain 
for valid working points of the technical system is determined. An efficient working point is 
determined using a Stratonovich equation and is used to design the technical system. 
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SUBSTITUTE SPECIFICATION 



METHOD AND ARRANGEMENT FOR DESIGNING A TECHNICAL SYSTEM 



[001] This application is the national phase under 35 U.S.C. § 371 of PCT International 
Application No. PCT/DE00/01035 which has an International filing date of April 4, 2000, 
which designated the United States of America, the entire contents of which are hereby 
incorporated by reference. 

Field of the Invention 

[002] The invention generally relates to a method and an arrangement for designing a 
technical system. 

Background of the Invention 

[003] In order to design a complex technical system, it is necessary to determine among a set 
of permissible working points (operating points, design parameters) at least one working 
point which permits effective running of the system. This frequently requires taking account 
of secondary conditions which limit the set of working points and therefore also need to be 
taken into account when searching for the working point. 

[004] In practice, the technical system is frequently designed with the aid of experience and 
the knowledge of one (or more) expert. It is disadvantageous in this case that the probability 
of error in the case of manual design rises when a certain complexity of the technical system 
is exceeded. 

[005] A method for optimization which follows an approach with the aid of which global 
optimization is carried out on the basis of a numerical solution of a stochastic differential 
equation is known from reference [1], listed at the end of the application and incorporated 
herein by reference. Secondary conditions are taken into account in this case as penalty 
functions. It is disadvantageous in this case that the range of permissible working points is 
left because 

a) information is thereby lost, 

b) permissibility is frequently more important than quality in a technical system, or 
it is not even possible to operate the system in an impermissible range, and 

c) penalty functions are difficult to prescribe or to dimension in practical use. 
[006] Reference [2] describes a method for local optimization under secondary conditions 
with the aid of projections, in the case of which nonlinear secondary conditions, in particular, 
are taken into account. 

[007] The concept of a Stratonovich equation (also: "Fisk-Stratonovich equation") is known 

1 
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from reference [3]. In this case, the Stratonovich equation denotes a type of stochastic 
differential equations with specific (see reference [3]) properties. 

SUMMARY OF THE INVENTION 
[008] An object of the invention is directed to permitting a design of a technical system 
which manages without separate penalty functions. 

[009] This object is achieved, for example, in accordance with the features of the 
independent patent claims. Developments of the invention follow from the dependent claims. 
[0010] In order to achieve the object, for example, a method for designing a technical system 
is specified in which the technical system comprises a target function which is influenced by a 
parameter vector having n parameters. Secondary conditions are prescribed with the aid of 
which a domain for valid working points of the technical system is determined. A projection 
operator is used to determine a Stratonovich equation which is solved numerically to 
determine an efficient working point. The determined efficient working point is used to 
design the technical system. 

[0011] As a development, the Stratonovich equation is determined in such a way that its 
solution varies in the domain of the valid working points. 

[0012] Another development is that the Stratonovich equation is solved numerically by way 
of a predictor-corrector method. A predictor can generally be described as a prediction 
variable which has, relative to the variable, to be predicted an error which is largely 
compensated by way of the corrector. 

[0013] In particular, the predictor can be determined with the aid of the trapezium rule. In 
this context, the trapezium rule for Riemann integrals is known, for example, from reference 
[4]. The trapezium rule for Stratonovich integrals follows from the definition of the 
Stratonovich integral itself. Thus, the Stratonovich integral is defined as a limiting value of 
the trapezium sums, and so the trapezium rule is a part of the definition of the Stratonovich 
integral. Trapezium sums are used to set up for the predictor an implicit equation which is 
solved approximately with the aid of a Newton method. 

[0014] One refinement resides in that step size control is additionally performed in the 
course of the global optimization. This takes place, in particular, byway of step size reduction 
if an excessively large approximation error is expected (value for the expected approximation 
error exceeds a prescribed bound), or by way of step size enlargement if a (very) small 
approximation error is expected (value for the expected approximation error is below a 
prescribed bound). 

[0015] It is also possible to carry out an additional step size reduction if an efficient working 
point determined using the above described method is used as a starting value for a 
subsequent local optimization method. It can be ensured thereby that an approximation with 
an increasingly smaller step size is performed in the vicinity of the local optimum, and the 

2 
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optimum is achieved with high accuracy thereby. 

[0016] Another refinement resides in that the Stratonovich equation has the following form, 
in particular: 

X e =X.-Jp(X.)-F(X.)-ds + Z (1) 

a 

with t > a and X a = x D 
where 

a denotes a starting instant (a > 0), 

x 0 denotes a permissible starting parameter vector, 

X t denotes a solution parameter vector at the instant t, 

P(X S ) denotes a projection matrix in the parameter vector X s , 

F(X S ) denotes an ascending direction in the parameter vector X s , 

Z denotes a random variable. 

[0017] Furthermore, the ascending direction can be determined by 

F= Vf(X s ), (2), 

where f denotes the target function. 

[0018] In an additional refinement, the random variable Z is determined by 
Z=E • Jp(X.)odB, (3), 

a 

where 

s denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian movement. 

[0019] It is known from reference [5] that a numerical solution of the stochastic differential 
equation 

t t 

X t = X a J Vf(X )ds + £ -JdBs.t^a (4) 
as 

a 

with an arbitrary starting parameter vector x 0 s R n (with R as the set of real numbers) and the 
starting condition X a = x G resulting therefrom results in a method for global optimization off 
in the parameter space. 



3 



New U.S. Application 
Docket No. 32860-000191 

[0020] If, in addition, secondary conditions are present in the form of equations and 
inequalities which limit the parameter space, the permissible parameter space is, in particular, 
initially reformulated, by way of introducing slack variables, such that it describes a manifold. 
In the above differential equation (4), it is only those components of the integrand which run 
tangential to the manifold that are taken into account. This is performed by multiplying by 
projection matrices. The result is the particular stochastic differential equation according to 
equation (1) (in particular seen in conjunction with equations (2) and (3), which is of the type 
of a Stratonovich differential equation. 

[0021] An additional development resides in that a secondary condition hi(x) < 0 which is 
present in the form of an inequality is transformed into an equation byway of a slack variable. 
This transformation has the form, in particular, of 

Hj_(x, s) = hi(x) + g(s), 



where 

g(s) denotes a function having a value range of all real values greater than or equal to 
zero (Rj ), 

s denotes the slack variable, and 
hi(x denotes a secondary condition hi(x) < 0. 

[0022] Furthermore, the function g(s) can be determined as: 
g(s) = s 2 (6). 

[0023] The solution of the equation (1) is a stochastic process 

{X t :t>a} (7) 
[0024] A sequence of permissible parameter vectors 

Xfc l /Xt 2 ' Xt 3 v " (a < ti < t 2 < t 3 ...) (8) 

is calculated which approximately represent points of a path of the stochastic process in 
accordance with the equation (7) at the instances t u t 2 , t 3 , .... An iteration step comes about in 
this case in the following way: let the parameter vectors x. ,x. , x, , ... ,x be already 

Z l t 2 fc 3 \ 

calculated, in which case x t then serves as starting point for the Stratonovich equation (1), 
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that is to say a: =tj. A parameter vector x t is calculated as a numerical approximation of 

i + 1 

a point on the path of the solution in accordance with equation (7), which already includes the 

points x fc ,x fc , x t , x t . The calculation is performed at the instant t i+ i, that is to say 

12 3 1 

t = tj+i, and the size of the iteration step is yielded as ti+i - tj. 

[0025] A subset of the sequence (8), in particular that with the smallest target function 
values, serves the purpose of determining the starting points for a subsequent local search for 
the optimum. 

[0026] The approximation of the Stratonovich equation (1) is performed for the (i+l)th 
iteration step using the following rule (with the aid of the trapezium rule); 

+ K P K) + P (*t 1+1 )) ■ (Et i+1 - B t J 

O) 

[0027] Equation (9) is an implicit equation for the random variable x^ . An 

i + 1 

approximation for a realization of the random variable X<_ can be determined 

i + 1 i + l 

therefrom with the aid of the Newton method. 

[0028] The point y t is impermissible in some circumstances. With the aid of a Newton 
l + l 

method, the corrector is now used to determine a point x^ in the vicinity of the point 

i + l 

y t . which is permissible and serves as starting value for a subsequent iteration step. 

[0029] Another refinement resides in that the design constitutes a new design, an adaptation 
or a control of the technical system. 

[0030] For a successful design, values which are efficient in each case are determined for the 
parameter vector and the working point. An efficient valuation of the parameter vector means 
that no parameter can be changed without thereby worsening a Q factor of the target function. 
[0031] A technical system can be a process engineering installation or other system which is 
to be designed or adjusted with reference to various parameters. In particular, the parameters 
of the parameter vector can be design parameters or operating parameters of the technical 
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system. Operating parameters characterize possible adjustable variables, whereas design 
parameters describe, in particular, physical dimensions of the technical system and can be 
adapted or varied during operation mostly only with a high outlay. 

[0032] The method described can optionally be used to redesign the technical system or to 
adapt an already existing technical system, hi both cases, this is a design (once as regeneration 
and once as adaptation) within the meaning of the present discussion. 

[0033] The technical system is realized or adjusted with the aid of the determined parameters 
within the scope of a further refinement. It is advantageous in this case that the parameters in 
a parameter vector which have been determined by the invention characterize a stable 
operating point, and that the adjustment of the system to this operating point ensures long- 
term reliable operation of the system/installation. 

[0034] To better achieve the object, moreover, an arrangement is specified for designing a 
technical system in the case of which a processor unit is provided which is set up in such a 
way that 

a) the technical system comprises a target function which can be influenced by a 
parameter vector having n parameters; 

b) secondary conditions are prescribed with the aid of which a domain for valid 
working points of the technical system can be determined; 

c) a Stratonovich equation is determined by means of a projection operator; 

d) the Stratonovich equation can be solved numerically and an efficient working 
point can be determined thereby; 

e) the efficient working point can be used to design the technical system. 

[0035] This arrangement is particularly suitable for carrying out the method according to the 
invention or one of its developments explained above. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0036] Exemplary embodiments of the invention are set forth and explained below with the 
aid of the drawing, in which: 

Figure 1 shows a sequence of operations of a method for designing a technical 
system; 

Figure 2 shows a sketch which illustrates a numerical solution of a Stratonovich 

differential equation; 
Figure 3 shows a sketch which illustrates a projection method; 
Figure 4 shows a processor unit. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0037] Figure 1 shows a block diagram which illustrates the sequence of operations of a 
method for designing a technical system. A target function and associated secondary 
conditions are provided in a block 101. The secondary conditions limit a space, which is 
determined by the number of the n parameters of the parameter vector, to a domain as part of 
this space. An efficient working point from this domain is determined as parameter vector by 
setting up and solving a Stratonovich equation by using projection operators which take 
account of the secondary conditions (compare blocks 102 and 103). Some of the working 
points found serve as starting points for a local search which, in turn, delivers an optimum 
point (efficient working point) (see block 104). In particular, the numerical solution to the 
Stratonovich equation includes a predictor-corrector method. 

[0038] Figure 2 illustrates the predictor-corrector method. Starting from a point x t which 



is located in the permissible working range L, a point is determined (application of 

l + l 

the Newton method to equation (9), for example) which is no longer situated in the 
permissible range L (predictor step 201). In order to transfer the point into the 

l + l 

permissible range again, a correction step 202 is determined which leads to the point 

x t. 
u i +1 

[0039] The interplay between predictor and corrector is repeated iteratively with a constant, 
decreasing or increasing step size. It may be noted in this case that the line 201 does not 
necessarily run tangentially. 

[0040] Figure 3 illustrates the projection method. Once again, the permissible range L is 
specified on which the point x t is located. The line 203 characterizes an ascending 

direction F(X t ), and a line 204 indicates a projected ascending direction P(X t ) • F(X t ). 
[0041] An example processor unit PRZE is illustrated in Figure 4. The processor unit PRZE 
includes a processor CPU, a memory MEM and an input/output interface IOS, which is used 
via an interface IFC in different ways: via a graphics interface, an output is visualized on a 
monitor MON and/or output on a printer PRT and/or output to any other type of output 
device. An input is performed via a mouse MAS and/or a keyboard TAST and/or any other 
type of input device. The processor unit PRZE also includes a data bus BUS which ensures 
the connection of a memory MEM, the processor CPU and the input/output interface IOS. 
Additional components, for example an additional memory, data memory (hard disk), 
scanner, etc., can also be connected to the data bus BUS. 

[0042] The invention being thus described, it will be obvious that the same may be varied in 
many ways. Such variations are not to be regarded as a departure from the spirit and scope of 
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the invention, and all such modifications as would be obvious to one skilled in the art are 
intended to be included within the scope of the following claims. 

[0043] List of references, the entire contents of each of which is hereby incorporated herein 
by reference: 

[1] K. Ritter, S. Schaffler: "A Stochastic Method for Constrained Global Optimization", 
SIAM 3. Optimization, Vol. 4, No. 4, pp. 894-904, 1994. 

[2] D.G. Luenberger: "Linear and Nonlinear Programming", Addison- Wesley Publishing 
Company, Massachusetts, pp. 334-337. 

[3] P. Protter: "Stochastic Integration and Differential Equations - A New Approach", 
Springer Verlag, New York, pp. 215-235. 

[4] J. Stoer: "Numerische Mathematik 1" ["Numerical Mathematics 1"], Springer- Verlag, 
Berlin, Heidelberg, 1994, pp. 138-144. 

[5] S. Schaffer: "Global Optimization Using Stochastic Integration"; Theorie und 
Forschung, Vol.340: Mathematik, Vol.3; Regensburg: Roderer, 1994; ISBN 3- 
89073-0; pp. 19-29. 
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Method and arrangement for designing a technical system 

THe Invention Relates to a method and an arrangement 
for designing a technical system. 

In order to design a complex technical system, it is 
necessary to determine among a set of permissible 
working points (operating points, design parameters)^ at 
least one working point which permits effective running 
of the system. This frequently requires taking account 
of secondary conditions which limit the set of working 
points and therefore also need to be taken into account 
when searching for the working point . 



In practice, the technical system is frequently 
designed with the aid of experience and the knowledge 
one (or more) expert. It is disadvantageous in this 
case that the probability of error in the case of 
manual design rises when a certain complexity of the 
technical system is exceeded. 



A method for optimization which follows an approach 

with the aid of which global optimization is carried 

out on the basis of a numerical solution of a a jl 

stochastic differential equation is known from^ ' Qa^ °^ 

Secondary conditions are taken into account in this i^^Pf 1 - 1 ^ 

case as penalty functions. It is disadvantageous in 

this case that the range of permissible working points t-eftrf'***- 

is left because o 

a) information is thereby lost, 

b) permissibility is frequently more important 
than quality in a technical system, or it is 
not even possible to operate the system in an 
impermissible range, and 

c) penalty functions are difficult to prescribe or 
to dimension in practical use. 
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^[2] describes a method for local optimization under 
secondary conditions with the aid of projections, in 
the case of which nonlinear secondary conditions, in 
particular, are taken into account. 

5 

The concept of a Stratonovich equation (also: "Fisk- 
Stratonovich equation") is known from) [3]. In this 
case, the Stratonovich equation denotes a type of 
stochastic differential equations with specific (seef^^ 
10 [3]) properties, _ \ 

The. — obg^et\ of the invention consists — ift permitting a 
design of a technical system which manages without 
separate penalty functions. 

This object is achieved ^in accordance with the features 
of the independent patent claims. Developments of the 
invention follow from the dependent claims. 

2 0 In order to achieve the object, j^a method for designing 
a technical system is specified in which the technical 
system comprises a target function which is influenced 
by a parameter vector having n parameters . Secondary 
conditions are prescribed with the aid of which a 

2 5 domain for valid working points of the technical system 
is determined. A projection operator is used to 
determine a Stratonovich equation which is solved 
numerically to determine an efficient working point. 
The determined efficient working point is used to 

30 design the technical system. 

As a development, the Stratonovich equation is 
determined in such a way that its solution varies in 
the domain of the valid working points. 



35 



Another development eonsi&fc-s- — Lu that the Stratonovich 
equation is solved numerically by mea^ls of a predictor- 
corrector method. A predictor can generally be 
described 
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as a prediction variable which has^ relative to the 
variable^ to be predicted an error which is largely 
compensated by meaftl of the corrector. 

In particular, the predictor can be determined with the 
aid of the trapezium rule. In this context, the 
trapezium rule for Riemann integrals is known, for 
example, from ; [4]. The trapezium rule for Stratonovich 
integrals follows from the definition of the 
Stratonovich integral itself. Thus, the Stratonovich 
integral is defined as a limiting value of the 
trapezium sums, and so the trapezium rule is a part of 
the definition of the Stratonovich integral. Trapezium 
sums are used to set up for the predictor an implicit 
equation which is solved approximately with the aid of 
a Newton method. 

One refinement consis^sj in that step size control is 
additionally performed in the course of the global 
optimization. This takes place, in particular, by wearfs 
of step size reduction if an excessively large 
approximation error is expected (value for the expected 
approximation error exceeds a prescribed bound) , or by 
mea«si of step size enlargement if a (very) small 
approximation error is expected (value for the expected 
approximation error is below a prescribed bound) . 

It is also possible to carry out an additional step 
size reduction if an efficient working point determined 
using the above described method is used as a starting 
value for a subsequent local optimization method. It 
can be ensured thereby that an approximation with an 
increasingly smaller step size is performed in the 
vicinity of the local optimum, and the optimum is 
achieved with high accuracy thereby. 

Another refinement cQnsigfcs in that the Stratonovich 
equation has the following form, in particular: 
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X t =X a -|p(X )-F(X s )-ds + Z (1) 

a 

with t > a and X a = x Q 
where 

5 a denotes a starting instant (a > 0) , 

x 0 denotes a permissible starting parameter 
vector, 

X t denotes a solution parameter vector at the 
instant t, 

10 P(X S ) denotes a projection matrix in the parameter 

vector X s/ 

F(X S ) denotes an ascending direction in the 

parameter vector X s , 
Z denotes a random variable. 



15 



20 



25 



Furthermore, the ascending direction can be determined 
by 

F = Vf(X B ) , (2) , 

where f denotes the target function. 

In an additional refinement, the random variable Z is 
determined by 



Z = £ 



Jp(X 3 )odB i (3) , 



where 

£ denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional 

3 0 Brownian movement. 

It is known from ^[5] that a numerical solution of the 
stochastic differential equation 

t t 

X t - X a J Vf (X )ds + £ -JdB s , t > a (4) 
a s 

a 
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with an arbitrary starting parameter vector x 0 e R n 
(with R as the set of real numbers) and the starting 
condition X a = x Q resulting therefrom con&fe-i-trutres a 
method for global optimization of f in the parameter 
5 space . 

If, in addition, secondary conditions are present in 
the form of equations and inequalities which limit the 
parameter space, the permissible parameter space is, in 

10 particular, initially reformulated, by mmt^ of 
introducing slack variables, such that it describes a 
manifold. In the above differential equation (4) , it is 
only those components of the integrand which run 
tangential to the manifold that are taken into account. 

15 This is performed by multiplying by projection 
matrices. The result is the particular stochastic 
differential equation according to equation (1) (in 
particular seen in conjunction with equations (2) and 
(3) , which is of the type of a Stratonovich 

20 differential equation. 

An additional development coasisfrs in that a secondary 
condition hi (x) < 0 which is present in the form of an 
inequality is transformed into an equation by means of 
25 a slack variable. This transformation has the form, in 
particular, of 

hi(x, s) = hi(x) + g(s), 

3 0 where 

g(s) denotes a function having a value range of 
all real values greater than or equal to zero 

(Rj), 

s denotes the slack variable, and 
35 hi(x denotes a secondary condition hi(x) < 0. 



Furthermore, the function g(s) can be determined as: 
= s 2 ( 6 ) . 
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The solution of the equation (1) constitofees a 
stochastic process 

{X t :t > a} (7) 

A sequence of permissible parameter vectors 

x*. ,x r /... (a < ti < t 2 < t 3 . . . ) (8) 

u l c 2 3 

is calculated which approximately represent points of a 

path of the stochastic process in accordance with the 

equation (7) at the instances ti, t 2 , t 3 , .... An 

iteration step comes about in this case in the 

following way: let the parameter vectors 

x ,x ,x ,...,x be already calculated, in which case 
t x t 2 t 3 t i 

x t then serves as starting point for the Stratonovich 
i 

equation (1), that is to say a: =ti. A parameter vector 

Xx. is calculated as a numerical approximation of a 

i + l 

point on the path of the solution in accordance with 
equation (7) , which already includes the points 

x. ,x , x , x. . The calculation is performed at the 
t l fc 2 t 3 i 

instant t i+ i, that is to say t = ti +x , and the size of 

the iteration step is yielded as ti +1 - ti. 

A subset of the sequence (8) , in particular that with 
the smallest target function values, serves the purpose 
of determining the starting points for a subsequent 
local search for the optimum. 

The approximation of the Stratonovich equation (1) is 
performed for the (i+l)th iteration step using the 
following rule (with the aid of the trapezium rule) : 
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X ti + 1 w X ti 

2 



K^ti) vf(x ti ) + p(x ti + 1 ).Vf(x t . +1 )) 
f( ? K) + P ( X t i + l))-K i+1 "BtJ 



(9) 



a 

m 

m 
a 

m 

s 

a 
in 

4= 
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Equation (9) is an implicit equation for the random 

variable x h . An approximation y<_ for a 

c i + 1 c i + 1 

realization of the random variable X h can be 

i + l 

determined therefrom with the aid of the Newton method. 

5 

The point is impermissible in some 

i + l 

circumstances. With the aid of a Newton method, the 

corrector is now used to determine a point x h in 

i + l 

the vicinity of the point y h which is permissible 

i + l 

10 and serves as starting value for a subsequent iteration 
step . 

resides 

Another refinement cQneir&fe-s in that the design 
constitutes a new design, an adaptation or a control of 
15 the technical system. 



For a successful design, values which are efficient in 
each case are determined for the parameter vector and 
the working point. An efficient valuation of the 
2 0 parameter vector means that no parameter can be changed 
without thereby worsening a Q factor of the target 
function. 



A technical system can be a process engineering 
25 installation or other system which is to be designed or 
adjusted with reference to various parameters. In 
particular, the parameters of the parameter vector can 
be design parameters or operating parameters of the 
technical system. Operating parameters characterize 
30 possible adjustable variables, whereas design 
parameters describe, in particular, physical dimensions 
of the technical system and can be adapted or varied 
during operation mostly only with a high outlay. 
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The method described can optionally be used to redesign 
the technical system or to adapt an already existing 
technical system. In both cases, 



a 

i» 

a 
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this is a design (once as regeneration and once as 
adaptation) within the meaning of the present discussion. 

The technical system is realized or adjusted with the 
aid of the determined parameters within the scope of a 
further refinement. It is advantageous in this case 
that the parameters in a parameter vector which have 
been determined by - me a ns — e€ the invention characterize 
a stable operating point, and that the adjustment of 
the system to this operating point ensures long-term 
reliable operation of the system/installation. 

in &^er to achieve the object, moreover, an 

arrangement is specified for designing a technical 
system in the case of which a processor unit is 
provided which is set up in such a way that 

a) the technical system comprises a target 
function which can be influenced by a parameter 
vector having n parameters; 

b) secondary conditions are prescribed with the 
aid of which a domain for valid working points 
of the technical system can be determined; 

c) a Stratonovich equation is determined by means 
of a projection operator; 

d) the Stratonovich equation can be solved 
numerically and an efficient working point can 
be determined thereby; 

e) the efficient working point can be used to 
design the technical system. 



This arrangement is particularly suitable for carrying 
out the method according to the invention or one of its 
developments explained above. 

Exemplary embodiments of the invention are set forth 
and explained below with the aid of the drawing, 
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in which: 



figure 1 shows a sequence of operations of a method 
for designing a technical system; 

5 

figure 2 shows a sketch which illustrates a numerical 
solution of a Stratonovich differential 
equation; 

figure 3 shows a sketch which illustrates a projection 
10 method; 



figure 4 shows a processor unit. 

Figure 1 shows a block diagram which illustrates the 

15 sequence of operations of a method for designing a 
technical system. A target function and associated 
secondary conditions are gi^en in a block 101. The 
secondary conditions limit a space, which is determined 
by the number of the n parameters of the parameter 

20 vector, to a domain as part of this space. An efficient 
working point from this domain is determined as 
parameter vector by setting up and solving a 
Stratonovich equation by using projection operators 
which take account of the secondary conditions (compare 

25 blocks 102 and 103) . Some of the working points found 
serve as starting points for a local search which, in 
turn, delivers an optimum point (efficient working 
point) (see block 104) . In particular, the numerical 
solution to the Stratonovich equation camp^isa^ a 

30 predictor-corrector method. 

Figure 2 illustrates the predictor-corrector method. 

Starting from a point x h which is located in the 

i 

permissible working range L, a point y- is 

u i + 1 
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determined ^application of the Newton method to 
equation (9)j) which is no longer situated in the 

permissible range L (predictor step 201) . In order to 

transfer the point y h into the permissible range 

u i + l 

again, 
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a correction step 202 is determined which leads to the 

point Xj. 

c i + 1 



The interplay between predictor and corrector is 
5 repeated iteratively with a constant, decreasing or 
increasing step size. It may be noted in this case that 
the line 201 does not necessarily run tangentially . 

Figure 3 illustrates the projection method. Once again, 
10 the permissible range L is specified on which the point 

x h is located. The line 203 characterizes an 

l + l 

ascending direction F(X t ), and a line 204 indicates a 
projected ascending direction P(X t ) • F(X t ). 

15 AJprocessor unit PRZE is illustrated in Figure 4. The 
processor unit PRZE c o mp r ise ^ a processor CPU, a memory 
£xf*E\ and an input /output interface IOS, which is used 
via an interface IFC in different ways: via a graphics 



20 



interface, an output is visualized on a monitor MON 
and/or output on a printer PRT/. An input is performed 
via a mouse MAS or a keyboard TASTVC The processor unit 
PRZE also -hars \ a data bus BUS which ensures the 
connection of a memory MEM, the processor CPU and the 
input/output interface IOS. Additional components, for 
25 example an additional memory, data memory (hard disk) 
scanner^ can also be connected to the data bus BUS. 



ft 
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e< -pzr — fcy i 
) in which the technical system comprise 

target function [which is/fT^Tfir^ a 
parameter -^ct^fh^^ c^fH^^^ 



^Jd) in whichj Secondary condi t ions^^^T^r escribed 

with the aid of whichj ^cSS^xn for valid 

working points of the technical system [is 

determined; ^ , it* 
(__c) m which) ^a Stratonovich equation [is determined 

by means of/jVj^oj ection operator; 
I d) in whichfjrae 4tratonovich equation fis solved"] 

numerically^ and'j an efficient working point Qs 

thereby de t e rmi n^ £=d z ^ 
Qe) in }fihi^^fie^Bft^%enfr working point (Is usedj to 

design the technical system. 

The method as claimed in claim 1, (jLn whichj the 
Stratonovich equation is determined (in) such (a wayO 
that its solution varies in (the) ^domain of (1:heQ 
valid^working points. 




The method as claimed in claim 1 [or s|, (Tn which) jJ^e*^ 
the Stratonovich equation is solved numerically by 
(means of/a predictor-corrector method. 

C^The method as claimed in claim 3, (in whic^jffhe 



predictor is determined with the aid of (thef^ 
trapezium rule. 




The method as claimed in claim 3 (^r (\n ^hich/Uf^M 
the corrector is determined by (means of thejpewton 



method. 
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The method as claimed in / one of 
claima,--rin whicHfstep size control is 
performed). 




the preceding 
s additionally 



m 
m 
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The method as claimed in claim 6, fin whictTl/the~ 
^a£>^tftep size control ^s performed in such a wayQ , ^ c. 

size reduction Q.S carried out^ if a 
sufficiently large approximation error is 
determined^ ^orj a,A^ 
b) tha^k :i:: ^tep size enlargement jis carried outfjif 
a sufficiently small approximation error is 
determined. 

The method as claimed in\( one of the preceding 
claims, B.n whichhThe Stratonovich equation has the 




following form: 



X t = X a - Jp(X s ) . F(X S ) - ds + Z^ [wit^Y 



t > a and X a = x Q 



r where"! n^A uJ^jc t n 
a denotes a starting instant (a > 0) , 
x 0 denotes a permissible starting parameter 
vector, 

X t denotes a solution parameter vector at the 
instant t, 

P(X S ) denotes a projection matrix in the parameter 
vector X s/ 

F(X S ) denotes an ascending direction, and 
Z denotes a random variable. 

The method as claimed in claimj[9|, (in whichjJlThe 
ascending direction is determined by; 

F = Vf(X fl ) , 

^^e^ f ^3enotingjj€he target function. 

10. The method as claimed in claim 8 (or ^, £in whichl/J^A 
the random variable is determined by 
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t 

Z = 8 • JP(X S ) o dB sir 

a 

8 denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian 
movement . 



J* 



The method as 
claimd^*""7Tri which] 



:Laimed in |one of the preceding 
the determined working point £is 
used^ as ^Starting point for a local optimization 
10 method. 

12 .L The method as clAimed in Tone of the preceding 

claim^^urT whicfijj a secondary condition which is 
present in the form of an inequality JJis 
15 transformed^into an equation ^by means ofJj'IfHslack 

variable s. 

13 V The method as claimed in claim 12, (in whic^Jthe 

secondary condition hi (x) , which is present in the 
20 form of the inequality hi(x)<0, is formulated as 

fii(x, s) = hi(x) + g(s) f 

C where] iS^m^ , , i 

g(s) denotes a function (having! fa value range of 
25 all real values greater than or equal to 

zero, 

s denotes the slack variable, and 

hi (x) denotes a secondary condition hi (x) < 0. 

30 14. The method as claimed in claim 13, ^in which] the 
function g(s) is determined as 



g(s) = s 2 . 



10 



15 



20 



GR 99 P 1730 



15 




J 



- 15 - 

J The method as claimed in (one of the preceding n 
claim^^^Tfn whicl^yxHe design /constitutes/ fit — new 
design, an adaptation ^orj^a control of the 
technical system . 

An arrangement for designing a technical system, C 6^H% ^ » 
QLti whic^ a processor unit (Is provided which is^ set 
up (in such a way th.foj& ^Wffr 



(ji)J the technical 
function wh: 



system [ comprisesj^pa target 



function which can be influenced by a parameter 
vector /having |pj JparametersA 1 a^M£~j 
rbi]feecondary conditions \ are prescribed with the 
aid of whicra^acfomain for valid working points 



G 



of the technical systemfcan be determined;/ 
c) {a — &t"ratOTiovich equation fis determined by means 



ofj^^pr^jection operator; 
^d^^T^e Stratonovich 
numerical 



equation 



be 



ng point ^can 



solved I 



eqi _ 
.ly and ^ an — efficient worki 
be determined thereby; o-wjcj 
[ e)mfhe j^grricient working point /can be usedj to 
design the technical system. 
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Abstract — — - - 



Method— and-a a? r ang ement_for__das igning— a- techn iea^^sys-trem 

A method for designing a technical system is specified 
in the case of which the technical system comp^rs^ a 
target function which is influenced by a parameter 
vector having n parameters. Furthermore, secondary 
conditions are prescribed with the aid of which a 
domain for valid working points of the technical system 
is determined. An efficient working point is determined 
by — mea*i« — % Stratonovich equation and is used to 
design the technical system. 
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Description 

| Method and arrangement for designing a teclinical system 



The invention relates to a method and an arrangement 
for designing a technical system. 

In order to design a complex technical system, it is 
necessary to determine among a set of permissible 
working points (operating points, design parameters) at 
least one working point which permits effective running 
of the system. This frequently requires taking account 
of secondary conditions which limit the set of working 
points and therefore also need to be taken into account 
when searching for the working point . 

In practice, the technical system is frequently 
designed with the aid of experience and the knowledge 
or one (or more) expert. It is disadvantageous in this 
case that the probability of error in the case of 
manual design rises when a certain complexity of the 
technical system is exceeded. 



A method for optimization which follows an approach 
with the aid of which global optimization is carried 
out on the basis of a numerical solution of a 
stochastic differential equation is known from [1] . 
Secondary conditions are taken into account in this 
case as penalty functions. It is disadvantageous in 
this case that the range of permissible working points 
is left because 

a) information is thereby lost, 

b) permissibility is frequently more important 
than quality in a technical system, or it is 
not even possible to operate the system in an 
impermissible range, and 

c) penalty functions are difficult to prescribe or 
to dimension in practical use. 
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[2] describes a method for local optimization under 
secondary conditions with the aid of projections, in 
the case of which nonlinear secondary conditions, in 
particular, are taken into account. 

The concept of a Stratonovich equation (also: u Fisk- 
Stratonovich equation") is known from [3] . In this 
case, the Stratonovich equation denotes a type of 
stochastic differential equations with specific (see 
[3] ) properties . 

The object of the invention consists in permitting a 
design of a technical system which manages without 
separate penalty functions. 

This object is achieved in accordance with the features 
of the independent patent claims. Developments of the 
invention follow from the dependent claims. 

In order to achieve the object, a method for designing 
a technical system is specified in which the technical 
system comprises a target function which is influenced 
by a parameter vector having n parameters. Secondary 
conditions are prescribed with the aid of which a 
domain for valid working points of the technical system 
is determined. A projection operator is used to 
determine a Stratonovich equation which is solved 
numerically to determine an efficient working point. 
The determined efficient working point is used to 
design the technical system. 

As a development, the Stratonovich equation is 
determined in such a way that its solution varies in 
the domain of the valid working points. 

Another development consists in that the Stratonovich 
equation is solved numerically by means of a predictor- 
corrector method. A predictor can generally be 
described 
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as a prediction variable which has relative to the 
variable to be predicted an error which is largely- 
compensated by means of the corrector. 

5 In particular, the predictor can be determined with the 
aid of the trapezium rule. In this context, the 
trapezium rule for Riemann integrals is known, for 
example, from [4] . The trapezium rule for Stratonovich 
integrals follows from the definition of the 

10 Stratonovich integral itself. Thus, the Stratonovich 
integral is defined as a limiting value of the 
trapezium sums, and so the trapezium rule is a part of 
the definition of the Stratonovich integral. Trapezium 
sums are used to set up for the predictor an implicit 

15 equation which is solved approximately with the aid of 
a Newton method. 

One refinement consists in that step size control is 
additionally performed in the course of the global 

20 optimization. This takes place, in particular, by means 
of step size reduction if an excessively large 
approximation error is expected (value for the expected 
approximation error exceeds a prescribed bound) , or by 
means of step size enlargement if a (very) small 

25 approximation error is expected (value for the expected 
approximation error is below a prescribed bound) . 

It is also possible to carry out an additional step 
size reduction if an efficient working point determined 

3 0 using the above described method is used as a starting 
value for a subsequent local optimization method. It 
can be ensured thereby that an approximation with an 
increasingly smaller step size is performed in the 
vicinity of the local optimum, and the optimum is 

35 achieved with high accuracy thereby. 

Another refinement consists in that the Stratonovich 
equation has the following form, in particular: 
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X t =X.- Jp(X,)-F(X s )-ds + Z (1) 

a 

with t > a and X a = x Q 
where 

5 a denotes a starting instant (a > 0) , 

x 0 denotes a permissible starting parameter 
vector, 

X t denotes a solution parameter vector at the 
instant t, 

10 P(X S ) denotes a projection matrix in the parameter 

vector X s , 

F(X S ) denotes an ascending direction in the 

parameter vector X s/ 
Z denotes a random variable. 



15 



20 



Furthermore, the ascending direction can be determined 
by 

F = Vf (Xs), (2)/ 
where f denotes the target function. 

In an additional refinement, the random variable Z is 
determined by 



25 



Z 



= E-Jp(X)odB, (3); 



where 

£ denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional 

3 0 Brownian movement. 

It is known from [5] that a numerical solution of the 
stochastic differential equation 
t t 

X t = X a j Vf (X )ds+£ -JdB s , t > a ( 4 ) 
a s 

a 
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with an arbitrary starting parameter vector x 0 s R n 
(with R as the set of real numbers) and the starting 
condition X a = x D resulting therefrom constitutes a 
method for global optimization of f in the parameter 
space . 

If, in addition, secondary conditions are present in 
the form of equations and inequalities which limit the 
parameter space, the permissible parameter space is, in 
particular, initially reformulated, by means of 
introducing slack variables, such that it describes a 
manifold. In the above differential equation (4) , it is 
only those components of the integrand which run 
tangential to the manifold that are taken into account. 
This is performed by multiplying by projection 
matrices. The result is the particular stochastic 
differential equation according to equation (1) (in 
particular seen in conjunction with equations (2) and 
(3), which is of the type of a Stratonovich 
differential equation. 

An additional development consists in that a secondary 
condition hi (x) < 0 which is present in the form of an 
inequality is transformed into an equation by means of 
a slack variable. This transformation has the form, in 
particular, of 

h ± (x, s) = hi(x) + g(s) , 

where 

g(s) denotes a function having a value range of 
all real values greater than or equal to zero 

(Rj), 

s denotes the slack variable, and 

hi (x denotes a secondary condition hi (x) < 0 . 



Furthermore, the function g(s) can be determined as: 

g(s) = s 2 (6) . 
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The solution of the equation (1) constitutes a 
stochastic process 

{X t :t > a} (7) 

A sequence of permissible parameter vectors 

x t ,Xj- f x f , ... (a < ti < t 2 < t 3 . . . ) (8) 

u l u 2 L 3 



10 is calculated which approximately represent points of a 
path of the stochastic process in accordance with the 
equation (7) at the instances ti, t 2 , t 3/ .... An 
iteration step comes about in this case in the 
following way: let the parameter vectors 

15 x ,x ,x ,,..,x, be already calculated, in which case 
u l u 2 u 3 r i 

x t then serves as starting point for the Stratonovich 
i 

equation (1) , that is to say a: =ti. A parameter vector 

x h is calculated as a numerical approximation of a 

i + l 

point on the path of the solution in accordance with 
2 0 equation (7) , which already includes the points 

x fc ,x fc / x i _ x t . The calculation is performed at the 
12 3 1 

instant ti +lf that is to say t = t± +1 , and the size of 
the iteration step is yielded as ti +1 - ti. 

25 A subset of the sequence (8) , in particular that with 
the smallest target function values, serves the purpose 
of determining the starting points for a subsequent 
local search for the optimum. 



30 



The approximation of the Stratonovich equation (1) is 
performed for the (i+l)th iteration step using the 
following rule (with the aid of the trapezium rule) ; 
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Vf k) + p(x tl+1 ).vf(x ti+ J 



(9) 



a 

si 
a 

fij 

P 

ru 
a 
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Equation (9) is an implicit equation for the random 
variable x fc . An approximation y^ for a 

l+l i+l 

realization of the random variable can be 

i + l 

determined therefrom with the aid of the Newton method. 
The point y t i + 1 is impermissible in some 

circumstances. With the aid of a Newton method, the 
corrector is now used to determine a point x h in 

i + l 

the vicinity of the point y t which is permissible 

i + l 

and serves as starting value for a subsequent iteration 
step. 

Another refinement consists in that the design 
constitutes a new design, an adaptation or a control of 
the technical system. 

For a successful design, values which are efficient in 
each case are determined for the parameter vector and 
the working point. An efficient valuation of the 
parameter vector means that no parameter can be changed 
without thereby worsening a Q factor of the target 
function. 

A technical system can be a process engineering 
installation or other system which is to be designed or 
adjusted with reference to various parameters. In 
particular, the parameters of the parameter vector can 
be design parameters or operating parameters of the 
technical system. Operating parameters characterize 
possible adjustable variables, whereas design 
parameters describe, in particular, physical dimensions 
of the technical system and can be adapted or varied 
during operation mostly only with a high outlay. 
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The method described can optionally be used to redesign 
the technical system or to adapt an already existing 
technical system. In both cases, 
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this is a design (once as regeneration and once as 
adaptation) within the meaning of the present discussion. 

The technical system is realized or adjusted with the 
5 aid of the determined parameters within the scope of a 
further refinement. It is advantageous in this case 
that the parameters in a parameter vector which have 
been determined by means of the invention characterize 
a stable operating point, and that the adjustment of 
10 the system to this operating point ensures long-term 
reliable operation of the system/installation. 



In order to achieve the object, moreover, an 
arrangement is specified for designing a technical 
15 system in the case of which a processor unit is 
provided which is set up in such a way that 

a) the technical system comprises a target 
function which can be influenced by a parameter 
vector having n parameters; 
2 0 b) secondary conditions are prescribed with the 

aid of which a domain for valid working points 
of the technical system can be determined; 
c) a Stratonovich equation is determined by means 
of a projection operator; 
25 d) the Stratonovich equation can be solved 

numerically and an efficient working point can 
be determined thereby; 
e) the efficient working point can be used to 
design the technical system. 



30 



This arrangement is particularly suitable for carrying 
out the method according to the invention or one of its 
developments explained above. 



35 Exemplary embodiments of the invention are set forth 
and explained below with the aid of the drawing, 
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in which: 

figure 1 shows a sequence of operations of a method 
for designing a technical system; 

5 

figure 2 shows a sketch which illustrates a numerical 
solution of a Stratonovich differential 
equation; 

figure 3 shows a sketch which illustrates a projection 
10 method; 

figure 4 shows a processor unit. 

Figure 1 shows a block diagram which illustrates the 
15 sequence of operations of a method for designing a 
technical system. A target function and associated 
secondary conditions are given in a block 101. The 
secondary conditions limit a space, which is determined 
by the number of the n parameters of the parameter 
20 vector, to a domain as part of this space. An efficient 
working point from this domain is determined as 
parameter vector by setting up and solving a 
Stratonovich equation by using projection operators 
which take account of the secondary conditions (compare 

2 5 blocks 102 and 103) . Some of the working points found 

serve as starting points for a local search which, in 
turn, delivers an optimum point (efficient working 
point) (see block 104) . In particular, the numerical 
solution to the Stratonovich equation comprises a 

3 0 predictor- corrector method. 

Figure 2 illustrates the predictor-corrector method. 

Starting from a point x^. which is located in the 

u i 

permissible working range L, a point y r is 

l + l 
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determined (application of the Newton method to 
equation (9) ) which is no longer situated in the 
permissible range L (predictor step 201) . In order to 

transfer the point into the permissible range 

i + 1 

again, 
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a correction step 2 02 is determined which leads to the 

point x^ 

L i + 1 

The interplay between predictor and corrector is 
5 repeated iteratively with a constant, decreasing or 
increasing step size. It may be noted in this case that 
the line 201 does not necessarily run tangentially . 

Figure 3 illustrates the projection method. Once again, 
10 the permissible range L is specified on which the point 
is located. The line 203 characterizes an 

ascending direction F(X t ), and a line 204 indicates a 
projected ascending direction P (X t ) * F (X t ) . 

15 A processor unit PRZE is illustrated in Figure 4. The 
processor unit PRZE comprises a processor CPU, a memory 
SPE and an input/output interface IOS, which is used 
via an interface IFC in different ways: via a graphics 
interface, an output is visualized on a monitor MON 

2 0 and/or output on a printer PRT. An input is performed 
via a mouse MAS or a keyboard TAST. The processor unit 
PRZE also has a data bus BUS which ensures the 
connection of a memory MEM, the processor CPU and the 
input/output interface IOS. Additional components, for 

25 example an additional memory, data memory (hard disk) 
or scanner can also be connected to the data bus BUS. 



10 
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A method for designing a technical system, 

a) in which the technical system comprises a 
target function which is influenced by a 
parameter vector having n parameters; 

b) in which secondary conditions are prescribed 
with the aid of which a domain for valid 
working points of the technical system is 
determined; 

c) in which a Stratonovich equation is determined 
by means of a projection operator; 

d) in which the Stratonovich equation is solved 
numerically and an efficient working point is 
thereby determined; 

e) in which the efficient working point is used to 
design the technical system. 

The method as claimed in claim 1, in which the 
Stratonovich equation is determined in such a way 
that its solution varies in the domain of the 
valid working points. 

The method as claimed in claim 1 or 2 , in which 
the Stratonovich equation is solved numerically by 
means of a predictor- corrector method. 

The method as claimed in claim 3, in which the 
predictor is determined with the aid of the 
trapezium rule. 

The method as claimed in claim 3 or 4, in which 
the corrector is determined by means of the Newton 
method. 
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The method as claimed in one of the preceding 
claims, in which step size control is additionally 
performed. 
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The method as claimed in claim 6, in which the 
step size control is performed in such a way 

a) that a step size reduction is carried out if a 
sufficiently large approximation error is 
determined, or 

b) that a step size enlargement is carried out if 
a sufficiently small approximation error is 
determined . 

The method as claimed in one of the preceding 
claims, in which the Stratonovich equation has the 
following form: 

t 

x t = X a - Jp(X s ) • F(X S ) • ds + Z with t > a and X a = x Q 

a 

where 

a denotes a starting instant (a > 0) , 
x 0 denotes a permissible starting parameter 
vector, 

X t denotes a solution parameter vector at the 
instant t, 

P(X S ) denotes a projection matrix in the parameter 
vector X s , 

F (X s ) denotes an ascending direction, and 
Z denotes a random variable. 

The method as claimed in claim 9, in which the 
ascending direction is determined by 

F = Vf (X B ) , 

f denoting the target function. 



The method as claimed in claim 8 or 9, in which 
the random variable is determined by 
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t 

Z = S • JP(X S ) o dB s 

a 

where 

s denotes a scaling factor, and 

{B s : s > a} denotes an n-dimensional Brownian 
5 movement . 



11. The method as claimed in one of the preceding 
claims, in which the determined working point is 
used as starting point for a local optimization 
10 method. 



12 . The method as claimed in one of the preceding 
claims, in which a secondary condition which is 
present in the form of an inequality is 
15 transformed into an equation by means of a slack 

variable s. 



13. The method as claimed in claim 12, in which the 
secondary condition hi (x) , which is present in the 
2 0 form of the inequality h ± (x)<0, is formulated as 



Si(Xf s) = hi(x) + g(s), 

where 

g(s) denotes a function having a value range of 
25 all real values greater than or equal to 

zero, 

s denotes the slack variable, and 

hi (x) denotes a secondary condition hi (x) < 0. 



30 14. The method as claimed in claim 13, in which the 
function g(s) is determined as 



g(s) = s 2 . 



The method as claimed in one of the preceding 
claims, in which the design constitutes a new 
design, an adaptation or a control of the 
technical system. 

An arrangement for designing a technical system, 
in which a processor unit is provided which is set 
up in such a way that 

a) the technical system comprises a target 
function which can be influenced by a parameter 
vector having n parameters; 

b) secondary conditions are prescribed with the 
aid of which a domain for valid working points 
of the technical system can be determined; 

c) a Stratonovich equation is determined by means 
of a projection operator; 

d) the Stratonovich equation can be solved 
numerically and an efficient working point can 
be determined thereby ; 

e) the efficient working point can be used to 
design the technical system. 
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German Language Declaration 


Prior foreign appplications 
PrloritSt beansprucht 




P^ritY Claimed, 


19919106.9 DE 


27.04.1999 


m □ 


(Number) (Country) 
(Nummer) (Land) 


(Day Month Year Filed) Yes No 
(Tag Monat Jahr etngereicht) Ja Nein 


(Number) (Country) 
(Nummer) (Land) 


□ □ 

(Day Month Year Filed) Yes No 
(Tag Monat Jahr elngereicht) Ja Nein 


(Number) "(Country) 
(Nummer) (Land) 


□ □ 

(Day Month Year Filed) Yes No 
(Tag Monat Jahr emgereicht) Ja Nein 


Ich beanspruche hierrnft gernass Absatz 35 der Zivil- 
prozessordnung der Vereinigten Staaten, Paragraph 
120, den Vorzug aller unten aufgefuhrten Anmel- 
dungen und falls derGegenstand aus jedem Anspruch 
dfeser AnmeSdung nicht in einer fruheren 
amerikanischen Patsntanmeldung laut dem ersten 
Paragraphen des Absatzes 35 der Zivilprozeliordnung 
der Vereinigten Staaten, Paragraph 122 offenbart ist, 
erkenne ich gemSss Absatz 37, Bundesgesetebuch, 
Paragraph 1.56(a) meine Pflicht zur Offenbarung von 
Informationen an, die zwischen dem Anmeldedatum 
der fruheren Anmeldung und dem nationalen Oder PCT 
international Anmeldedatum dieser Anmeldung 
bekannt geworden sind. 


1 hereby claim the benefit under Title 35. United States 
Coda §120 of any United States application(s) listed 
below and, insofar as the subject matter of each of the 
claims of this application is not disclosed in the prior 
United States application in the manner provided by 
the first paragraph of Title 35, United States Code, 
§122, I acknowledge the duty to disdose material 
information as defined in TrtJe 37, Code of Federal 
Regulations, §1 .56(a) which occured between the filing 
date of the prior application and the national or PCT 
international filing date of this application. 


PCT/DEOO/01035 


04.04.2000 


anhSnqig sendinq 


(Application Serial No.) 
(An meldesa riennummer) 


(FiJing Data D, M.Y) 
(Anmeldedatum T. M, J) 


(Slntus) (Status) 
(patonliert, anhanglg, (patented, pending, 
aufgegeben) a ban dan ed) 


(Application Senaf No,) 
CAnmeldBsenennummQrj 


(Filing Date D,M.Y) 
(Anmeld&daltlm T, M. J) 


(Statue) (Status) 
(pGtentjert, anhangig, (patented, ponding, 
aufgcbsn) abandoned) 


Ich erklaYe hiermrt, dass a He von mir in der vorliegen- 
den Erklarung gemachten Angaben nach meinem 
besten Wissen und Gewissen der voilen Wahrhert 
entsprechen, und dass ich diese eidesstattiiche Erkla- 
rung in Kenntnis dessen abgebe, dass wissentl ich und 
vorsatziich falsche Angaben gem&ss Paragraph 1001, 
Absatz 18 der Ziviiprozessordnung der Vereinigten 
Staaten von Amerika mit Geldstrafe belegt und/oder 
Gefangnts bestraft werden koennen, und dass derartig 
wissentlich und vorssitzlich falsche Angaben die Gul- 
tigkeit der vorliegenden Patentanmeldung Oder ernes 
darauf erteitten Patentes gef&hrden konnen. 


I hereby declare that a!f statements made herein of my 
own knowledge are true and that all statements made 
on information and belief are believed to be true, and 
further that these statements were made w/th the 
knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, 
under Section 1001 of Title 18 of the United States 
Code and that such willful false statements may 
jeopardize the validity of the application or any patent 
issued thereon. 
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German Languagi 

VERTRETUNGSVOLLMACHT: Als benannter Erfinder 
beauftrage ich hiermrt den nachstehend benannten 
Patentanwalt (Oder die nachstehend benannten 
Patentanwaite) und/oder Patent-Agenten mit der 
Verfo!gung der vorllegenden Patentanmeldung sowie 
mit der Abwicklung alier damrt verbundenen Geschafte 
vor dem Patent- und Warenzeichenamt: (Nsme und 
Rogistrationsnummer anfQhren) 



Declaration 

POWER OF ATTORNEY: As a named inventor, I 
hereby appoint the following attorney(s) and/or 
agent(s) to prosecute this application and transact afl 
business In the Patent and Trademark Office 
connected therewith (list name and registration 
number) 



And 1 hor«by appoint 

Customer No. 02292 



Telefongesprache bitte richten an: Direct Telephone Calls to: (name and telephone 

(Name und Telefonnummer) number) 



Ext. 



Postanschrift: Send Correspondence to: 

Birch, Stewart, Kolasch & Birch, LLP 
8110 Gatehouse Road / Suite 500 East 22042 Falls Church, Virginia 
Telephone: +1 703 205 8000 and Facsimile +1 703 205 8050 

or 

Customer No. 02292 



VollarName dos einzigen Oder ureprOnghchen Erfinders: 

ANNEUE STOEHR _ 


Full name of soJe or first inventor 

ANNELIE STOEHR 


Untarschrift des Erflnders Datum 

A- — vSfoU /R/io.d 


Inventor's sig nature c%\q 

A~<JLz. Sfc'U Al./o.OA 


Wqhnsrfe 

MUENCHEN, DEUTSCHLAND ^>^~ 


Residence 

MUENCHEN, GERMANY 


Staatsang anting keit 

DE 


CHlzenship 

DE 


HEIDELBERGER STR.12 


Post Office Addess 

HEIDELBERGER STR.12 


80804 MUENCHEN 


80804 MUENCHEN 




Full name of second join! inventor, rf any: 


Unterechnrl des Erfinders Datum 


Second Inventor's stgnalure Date 




Residence 


Slaatsanc(*h6rigke/t 


Cfenship 


PoatansohfrR 


Post Office Address 







Fafle von dntten undweiieren Miterfindern sngebon). 



subsequent joint inventors). 
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